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Abstract. We explore the imprint of the cosmological hydrogen recombination lines on the power spectrum of the cosmic 
microwave background (CMB). In particular, we focus on the three strongest lines for the Balmer, Paschen and Brackett series 
of hydrogen. We expect changes in the angular power spectrum due to these lines of about 0.3/iK for the Ho- line, being maxi- 
mum at small angular scales {( ~ 870). The morphology of the signal is very rich. It leads to relatively narrow spectral features 
(Av/v ~ 10 _1 ), with several regions in the power spectrum showing a characteristic change of sign of the effect as we probe 
different redshifts or different multipoles by measuring the power spectrum at different frequencies. It also has a very peculiar 
dependence on the multipole scale, connected with the details of the transfer function at the epoch of scattering (see movies 
illustrating the effect at http://www.mpa-garching.mpg.de/~jalberto l. In order to compute the optical depths for these transi- 
tions, we have numerically evolved the populations of the levels of the hydrogen atom during recombination, simultaneously 
treating the evolution of helium. For the hydrogen atom, we follow each angular momentum state separately, up to the level 
n = 10. Foregrounds and other frequency dependent contaminants such as Rayleigh scattering may be a important limitation 
for these measurements, although the peculiar frequency and angular dependences of the effect that we are discussing might 
permit us to separate it. Detection of this signal using future narrow-band spectral observations can give information about the 
details of how the cosmic recombination proceeds, and how Silk damping operates during recombination. 
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1. Introduction 

Existing and planned experiments devoted to the measurements 
of Cosmic Microwave Background (CMB) angular fluctua- 
tions will reach an unprecedent high sensitivity of measure- 
ments. Planck, SPT, ACT and QUEST will achieve sensitivi- 
ties several orders of magnitude higher than the amplitude of 
the acoust ic peaks which were the dream of theorists in the 
seventies JSunvaev & Zeldovichl Il970t IPeebles & Yul Il970t 
lDoroshke^cliet^lTTl978T) '. and are observed with good preci- 
sion now by Boomerang, Maxima, DASI, VS A, CBI, and more 
recen tly with very high accuracy by WMAP (Bennett et al., 
l2003l) . 

The experimental progress demonstrates that we are enter- 
ing the era of "precision cosmology", and many of the effects 
that were obvious for theoreticians but of small amplitude be- 
come now accessible to experimentalists. 
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One of such problems is the direct detection of the hydro- 
gen or helium lines from the epoch of recombination of hy- 
drogen in the Universe at redshift z ~ 1000. In Russia, all 
the study of cosmological recombination grew up from the 
question of Vladimir Kurt: "where are the Ly-o- line photons 
from reco mbinatio n in the Univers e? ". This question forced 
IZeldovich. Kurt & Svunvaevl ( ll968l hereafter, ZKS68) to study 
the detailed physics of recombination, and to show that two- 
photon decay of 2s level of hydrogen is more important than 
the escape of redshifted Ly-a photons due to expansion of 
the Universe in determining the rate at which recombination 
occurs. It was shown that there are strong distortions of the 
CMB spectrum due to two-photon decay, but unfortunately 
these distortions are practically unobservable because they lie 
in the distant Wien region of the CMB spectrum. Peebles 
dl968l hereafter P68), and ISeager. Sasselov & Scottl dl999l 
2000, he reafter SSS99 and SSS00) came to the same conclu- 
sion. IDubrovichI i 19751) proposed to look for CMB spectral 
distortions due to transitions between higher levels (Balmer, 
Paschen and higher series), which should be broadened due to 
the delayed recombination, and many papers were devoted to 
the computation of these very small but extremely interesting 
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effects ( Grachev & Dubro vichlll99ll IdelT Antonio & RvbickH 
ll993UDubrovich & Shakhvorostoval I2004 . Although several 
computations of the int ensity of these lines gave co ntra- 
dictory results (see e.g. [ d ell 7 Anto nio & Rvbickil Jl993l) and 
iDubrovich & Shakhvorostovari2004l) '). it is obvious that the 
corresponding distortions of the CMB spectrum are extremely 
weak, and that AI/B y does not exceeds the level of 10~ 7 , which 
is two orders of magnitude below the sensitivity of COBE- 
FIRAS experiment, but certainly will be the goal of future more 
precise experiments. 

Another physical process with information about the epoch 
of recombination is Rayleig h scattering on neutra l hydrogen 
and helium atoms JYu. Soergel & Ostrikerl l200lh . This pro- 
cess, due to the characteristic dependence of the scattering 
cross-section on wavelength (cr R oc y 4 ), has a very strong fre- 
quency dependence of the power spectrum, C[. Unfortunately, 
this frequency dependence is rather similar to the spectrum 
of dust emission from local dust and bright extragalactic star- 
forming galaxies, and the amplitude of the effect becomes im- 
portant in frequencies where the dust emission is significant 
(v ^ 200 GHz). This will make a rather difficult task to sepa- 
rate it from nearby and bright foreground dust contribution. 

In this paper, we concentrate on the detection of hydro- 
gen lines, and we propose a different method to detect the 
consequences of the overpopulated levels of hydrogen atom 
at the epoch of recombination. We propose to look for fre- 
quency dependent angular fluctuations of the CMB in differ- 
ent angular scales and at frequencies in the vicinity of the red- 
shifted Balmer (and higher series) lines. Here, we estimate how 
the small optical depths connected with these transitions in- 
fluences the angular distribution of the CMB, arising mainly 
due to Thomson scattering on free electrons. We show that one 
should expect changes in the angular power spectrum due to 
these lines of about 0.3 fiK, being maximum at small angu- 
lar scales ({ w 870). These numbers are at the level of the 
planned sensitivity of the Planck satellite but for large angu- 
lar scales, and certainly should be achieved by future ground- 
based instruments looking for polarization of the CMB, and 
future space missions like CMBPol. It is important to mention 
that present day technology permits one to achieve sensiti vities 
at least 50 times better than that of Planck JChurchlEooH) . The 
maximum signal that we are discussing corresponds to the red- 
shifted Ha recombination line, which unfortunately is located 
in a frequency region where the contribution from dust emis- 
sion is very important, and where the Rayleigh scattering is 100 
times larger in Q (10 times in temperature). 

Coherent scattering in the Ho- line produces a AT /To w 
10" 7 signal in the angular distribution of the CMB, whereas 
spectral distortions in the same line are close to 10~ 8 - 10~ 7 . 
Currently, experimentalists put their main effort in observing 
angular fluctuations. Therefore, the signal under discussion in 
this paper might be observed sooner or later. 

This effect is small, but it leads to relatively narrow spec- 
tral features in the power spectrum (Av/v ~ 10 _1 ). In addi- 
tion, there are several regions in the angular power spectrum at 
high multipoles where the sign of <5Q (=Cj, lne - Q) connected 
with this process depends on the wavelength, i.e. it changes 
from positive to negative when the frequency is shifted from 
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Fig. 1. Distortion of the power spectrum due to coherent scat- 
tering in the Ha line, as a function of the observing frequency 
today. The maximum distortion (in absolute value) is obtained 
at v = 452 GHz for £ = 873. We show that shifting the observ- 
ing multipole in a region of At = ±70 around the maximum, 
we still have the signal within the 70% of the peak intensity. 

350 GHz to 450 GHz (Fig.UJ. The effect also has a very pe- 
culiar dependence on I at all frequencies (Fig. and these 
peculiarities in frequency and angular scale dependence might 
permit future observers to separate coherent scattering on Ha 
during the epoch of recombination from foregrounds and even 
from the cosmic variance of the effect of Rayleigh scattering. 

The discovery of these spectral features in the CMB angular 
distribution will permit us to measure directly the tail of the re- 
combination process and its position with good precision, and 
may give additional information about cosmological parame- 
ters of the Universe and their evolution with time. 

This may also give direct spectroscopic evidence of recom- 
bination. If Rayleigh scattering probes the existence of neutral 
hydrogen immediately after the peak of the visibility function, 
the profile of the Balmer lines features in the power spectrum 
gives us the possibility to check the dynamics of recombination 
when the electron density changes from n e /nu a; 1 to 10~ 2 , 
and correspondingly, the optical depth in Her and the popula- 
tion of the second level are traced. The characteristic physi- 
cal length producing fluctuations in Ha in that epoch corre- 
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, 2. Distortion of the power spectrum due to scattering in the 
line as a function of angular scale. 



eration of new angular fluctuations is connected with just re- 
distribution of radiation over angles, so it has no relation to the 
difference between excitation temperature and color tempera- 
ture of radiation. Optical depths in all lines, including Balmer 
lines, is very small (e.g. t ~ 5 X 10~ 5 for Ha line at z = 10 00). 
But as it was shown by iHernandez-M onteag udo & Sunvaevl 
d2005l hereafter HMS05) and BHMS04, correlation effects 
with the existing radiation fluctuation field amplifies the ef- 
fect in such a way that 6Ccs are effectively proportional to r v , 
instead of tJ. It should be noted that the signal that we are 
discussing is smaller than the cosmic variance level, but as is 
explained by HMS05 and BHMS04, multi-frequency observa- 
tions permit one to avoid the limit imposed by the cosmic vari- 
ance associated with the intrinsic CMB fluctuations, and reveal 
signals below this threshold. 

Finally, it is important to mention that these detections 
strongly rely on a accurate cross-calibration between different 
channels of a given experime nt. However, there is a p ossibility 
to use the thermal SZ effect dSunvaev & Zeldovichl 1 19721) on 
clusters of galaxies and sum of blackbodies appr oach for such a 
calib ration. Both methods are discussed by Chluba & Sunvaevl 
d2004l) . 

The outline of the paper is as follows. In Section 2, we re- 
view the basic equations which describe the cosmic recombina- 
tion process, and we present a code that we have developed to 
compute the optical depths in the lines for the required tran- 
sitions. In Section 3, we show how to derive the 5C( quan- 
tities for a given line, and present our results for the consid- 
ered transitions. Section 4 presents a discussion about the fore- 
ground contamination and the amplitude of the effect compared 
to Rayleigh scattering. Conclusions are given in Section 5. 



sponds to smaller scales than those for Rayleigh scattering. It 
is also important that for the hydrogen lines, each observed fre- 
quency corresponds to a given redshift z, but the 'picture' from 
Rayleigh scattering is the same for all frequencies (there is only 
a change in amplitude). Thus, when the line is placed inside re- 
combination, the shape of the angular power spectrum for these 
two effects differs strongly. 

Computations of the 5Cc connected with recombinational 
lines are direct, because detailed estimates of the recombina- 
tion process show (SSS99) that the assumption of ZKS68 about 
the Saha equilibrium of all levels above level two is valid with 
good precision. This gives us the possibility to find the opti- 
cal depth in each line of interest, and include it as an additional 

■ ■ 1 

opacity coefficient in the CMBFast code ( Seliak & Zaldarriaga, 

Il996l) . as it was done for the lines i n the case of reso- 
nant scatter ing by IZaldarriaea & Loebl d2002l) (lithium dou- 
blet), and dBasu. H ernandez-Monteagud o & Sunvaevl ^004, 
hereafter BHMS04) for fine structure lines of neutral atoms and 
ions of carbon, nitrogen and oxygen. Let us mention here that 
the computations of the line intensities (i.e. spectral distortions 
of CMB) require the detailed calculation of tiny deviations 
from Saha equilibrium, which are connected to the process of 
recombination controlled by 2s-level decay. For our purposes, 
it is enough to use the equilibrium distribution because the gen- 



2. Cosmological recombination process 

In this section, we review the basic equations describing the 
process of cosmic recombination of hydrogen and helium, and 
we show how to use them to infer the optical depths for the 
recombination lines. Since the first computations by ZKS68 
and P68, many refinements have been introduced (for a re- 
view, see SSS00). SSS00 has the most detailed calculations 
to date, including 300 levels for the hydrogen atom, 200 for 
Hel and 100 for Hell. When compared with the standard "ef- 
fective three-level computation", their calculation results in a 
roughly 10% change in the electron fraction at low redshift for 
most of the cosmological models, plus a delay in Hel recom- 
bination. All these results were incorporated in the RECFAST 
code (SSS99), which consists in a modification of the effective 
three-level model to reproduce the new values. 

For this paper, we are interested in computing the optical 
depths associated with the hydrogen recombination lines of the 
Balmer, Paschen and Brackett series. We will not di scuss here 
the helium lines ("see iDubrovich & Stoivarovl (ll997) for a dis- 
cussion on the spectrum of primordial Hel and Hell recombi- 
nation lines). In order to obtain the optical depths, we need to 
know the populations for all levels in the hydrogen atom which 
contribute to the transitions we want to study. We will do this 
following two different approaches. 
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On the one hand (method I), we will follow SSSOO and we 
will evolve the level populations for a hydrogen atom with 10 
levels, separating levels both by principal quantum number n 
and angular momentum /. In method II, we will adopt the ion- 
isation fraction from RECFAST and compute, using the Saha 
equation relative to the continuum, the popula tion of the ex- 
cited levels (e.g. iLiubarskii & Sunvaevl dl983l) ). This method 
is a valid approximation for those levels above the second one 
and for redshifts z S 900, but as we will see, this is indeed our 
range of interest because at lower redshifts the visibility func- 
tion is very small. 

2.1. Method I: Detailed follow-up of level population 

The detailed formalism and the equations describing the evo- 
lution of the population of the levels of hydrogen and helium 
atoms with cosmic time is presented in SSSOO. Here, we enu- 
merate the assumptions of our work and the differences to that 
paper. 

a) We follow in detail all sub-levels of the hydrogen atom, 
taking into account the quantum number / for angular 
momentum, and the principal quantum number n up to 
a maximum value of n max = 10. Thus, each state is 
labeled with two numbers, (n, I), and we have a total 
number of /V = n max (n max + l)/2 = 55 bound levels. 
It should be noted that SSSOO do not separate different 
levels according to the quantum number I, but thi s is in- 
cluded in the computation of dell ' Antonio & Rvbicki 
GUI- 

We decided to follow all separate I levels to 
cross-check our r esults for the intensity of s pectral lines 
with those from IdelF Antonio & R ybickil Jl993l) and 
Dubrovich & Shakhvorostova (2004). 

b) Only radiative rates will be taken into account when 
solving the set of equations, because collisional pro- 
cesses are not important in the early Universe (see 
ZKS68, P68 and SSSOO). 

c) For helium, we will not make a detailed follow-up of 
all levels, and we only evolve the equation for its total 
ionization fraction (see RECFAST code, SSS99). 

d) It will be assumed that the radiation field is a perfect 
blackbody (COBE/FIRAS showed that deviations from 
blackbody are sma ller than 10~ 4 JMather et ail 1 1994 
iFixsen etall IT996)). so the equation for the evolution 
of the radiation field will be omitted. This approxima- 
tion is not valid in the case of the Lyman series for hy- 
drogen, where the spectral distortions are strong, but it 
is good for Balmer series and above. 

e) We will use th e Sobolev escape probability method 
llSobolevH 19601) to deal with the evolution of resonance 
lines. This will decouple the evolution of the radiation 
field from the evolution of the level population. We note 
that this method was developed for a completely differ- 
ent problem. For a derivatio n of this method in a co s- 
mological situation, see lDubrovich & Grachev! (12004). 

Values for the physical constants for hydrogen and he- 
lium were taken from SSS99 and SSSOO. The values for 



the oscillator strengths and the corresponding Einstein coef- 
ficients for ins tantaneous emission were computed following 
iGreen. Rush & Chandler! Jl957l) . The value for the two photon 
decay transition A2 S ,i s is essential, given that this is the dom- 
inant mechanism driving the cosmological recombinati on. W e 
adopt here the latest value, 8.22458s -1 , from lGoldmar] Jl989). 

The photoionization/recombination rates were computed in 
two different ways. The f irst one was by using the photoioniza- 
tion cross-sections from Storey & Humme J dl99lh : iHummerl 



(1994). The second one uses the analytical expressions from 
Burgess ( 1958), which are valid for small values of the energy 
of the ejected electron. In both cases, radiative recombination 
rates included spontaneous and induced transitions. Our final 
results were obtained using the second method, which is com- 
putationally faster than the first one. Nevertheless, we checked 
that both methods give similar results for the populations of the 
levels. 

Once the populations of all levels are obtained at all red- 
shifts, the optical depth associated with any permitted transi- 
tion j — > i connecting the levels i - (n, I) and j = (m, I') (with 
m > n) can be computed using the Sobolev formula as 



Aji^j[ni(gi/gj)-nj] 
SttH(z) 



(1) 



where Ajj is the Einstein coefficient for instantaneous emission, 
gi is the degeneracy of the energy level, and H(z) is the (time 
dependent) Hubble constant. The total optical depth r„ m asso- 
ciated with a particular line m — > n is easily derived by adding 
the contribution from all the permitted transitions connecting 
the sublevels with those principal quantum numbers. 

In Fig.EJwe show our result for the electron density evolu- 
tion (x e = n e jnu) as a function of redshift compared with the 
standard calculation using RECFAST. These results are consis- 
tent with those obtained by SSSOO when they use n — 10 levels, 
but it should be noted that in our calculation we do a detailed 
follow-up of the population of the different angular momen- 
tum states. When more levels are included in the computation, 
SSSOO show that the residual electron density at low redshifts 
becomes smaller, and for n ^ 50 it converges as the atom be- 
comes complete in terms of energy levels. However, for the 
purposes of this paper, it is enough to consider a 10-level hy- 
drogen atom to achieve good precision in the region of interest. 

2.2. Method II: Saha equation from the continuum 

As suggested by ZKS68, and confirmed with the computations 
of SSSOO, the approximation of equilibrium of all levels above 
the second one yields very good results, and deviations greater 
than 10% are only found at redshifts z 800. Thus, we can 
derive (to a good approximation) the population of the level in 
two steps. We first use the standard RECFAST computation, 
and derive the evolution of the electron density. From here, and 
using the Saha equation relative to the continuum, we can de- 
rive the population of all levels above the second one, and can 



1 Note that in our notation, n H = «o(l + z) is the total number 
density of hydrogen, so n H = n H1 + n p . 
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Fig. 3. Multilevel hydrogen recombination as computed in this 
paper (using 10 levels for quantum number n, and follow- 
ing independently all the states of angular momentum), com- 
pared with the standard 'effective three level' computation from 
RECFAST. It is shown x e (= n e lnu) as a function of redshift in 
the vicinity of the epoch of cosmic recombination for the cos- 
mological m£dej_wjthparameters - 0.044, Q. to , = 1 and 
h = 0.71 JBennett et alll2003l) . 



Fig. 4. Optical depths for several lines considered in this pa- 
per, as a function of redshift z. All values were computed using 
method I. For the Ha line, it is also shown the computation 
using the Saha approximation (method II), to show that with 
this method, deviations with respect to the detailed calculation 
(method I) only appear when the Universe is practically trans- 
parent, which corresponds to the high frequency wing of the 
line and where the expected effect is not important. 



compute all the optical depths for B aimer lines and higher se- 
ries. It should be noted that this approximation fails for red- 
shifts below z ~ 900 because most of the atoms have recom- 
bined and the transition rates are not enough to keep equilib- 
rium between high levels. 

Thus, the population of the level n is derived in this second 
method as 

(2mn e k B T e yi L 2 

where the state i = (n, I) has excitation energy below the con- 
tinuum^, = 13.598/n 2 eV, and g,- is the statistical weight of 
the state. From that expression, we can derive the absorption 
coefficient for a transition / — > i, and from here and using the 
Sobolev approximation, the optical depth for each line at red- 
shift z (equation^. 

2.3. Optical depth in the hydrogen recombination lines 

Although we have computed many more transitions, in this pa- 
per we shall concentrate on those lines whose expected contri- 
bution is largest. Thus, we investigate here the three first tran- 
sitions (labelled as m — > n) for the Balmer (n = 2), Paschen 
(n = 3) and Brackett (n = 4) series. Table ^presents the wave- 
lengths and oscillator strengths for these lines. For illustration, 
we show in the last two columns the redshifted frequency and 
the optical depth in the line evaluated at redshift z = 1000, and 
for a cosmological model with parameter values taken from 
iBennett et ail i2003l) . i.e. Q, , = 1, Q. h = 0.044 and reduced 
Hubble constant h =0.71. Throughout this paper we have used 
these values, unless otherwise stated. 

Figure 0] shows the optical depth in the lines (j nm ) derived 
with our method I, and for this cosmological model. Using the 



method II, we obtain these same values for the optical depths, 
but only for redshifts above z ~ 900. To illustrate this point, 
we also display in Fig.|4]the optical depth for the Ha line us- 
ing the method II (dotted line). Method II reproduces the non- 
equilibrium computation above z ~ 900, but fails below this 
point, showing a divergent behaviour. The reason for this was 
pointed out above. At these redshifts, the radiation field has not 
enough photons to maintain the equilibrium of the second level 
with the continuum. The photoionization rates become smaller 
than the photorecombination rates, and the Saha equilibrium 
between the second level and the continuum is no longer valid 
(e.g. ZKS68). Moreover, the populations of 2s and 2p levels 
start to show strong departures from their equilibrium ratio. 
Given that the cosmological recombination is proceeding much 
more slowly than expected from a Saha recombination, method 
II (which uses the Saha equation relative to a continuum level 
which has been computed using the effective three level cal- 
culation) predicts very high values for the population of the 
second level in this redshift range. This is the reason for the 
divergent behaviour at z ^ 900 of the Ha optical depth com- 
puted with method II. Although we have used the exact (non- 
equilibrium) computation for all the results presented in this 
paper, we wanted to show that using a very simple approxima- 
tion (method II) it is possible to correctly estimate the values 
for the optical depths in the range of interest (as we show below, 
the peak of the effect we are investigating occurs at z ~ 1010). 

In Fig. |4] we finally note that for Balmer transitions, the 
shape of the optical depth as a function of frequency is similar 
for all lines. This can be easily understood because the popu- 
lation of the second level is much larger than the population of 
higher levels, and thus the optical depth is directly proportional 
to this population {n% s + ri2 P ) for all Balmer lines. 
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Table 1. List of the strongest hydrogen recombination lines studied in this paper. We present the three strongest lines for the 
first three series of hydrogen atom between excited states: Balmer (m — > 2), Paschen (m — > 3) and Brackett (m — > 4). For each 
transi tion m — > h, we show its name, the rest frame wavelength, and the oscillator strength (these values have been taken from 
iGreen. Rush & Chandlerldl957l) '). For illustration, we also present in the fifth column the redshifted frequency observed today if 
emission took place at z = 1000. The last column shows the optical depth in the line at that redshift (z - 1000), as computed 
using our code. 
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Fig. 5. Dependence of the optical depth with the redshift for the 
Ha line. For comparison, we also show the Thomson visibility 
function CWC^) = (dTj/dr]) exp(-T7-)), the function exp(-Tr), 
and the function TH a exp{—Tj). To display all of them on the 
same scale, they were rescaled by the indicated factors. 



3. Imprint of hydrogen recombination lines on the 
CMB 

The imprint o f line transitions on the C MB spectrum has been 
examined by Zaldarriag a & Loebl(l2002l) for the case of lithium 
recombination, and for other ions and atoms in BHMS04. We 
extend these works here and we consider the case in which the 
coherent scattering occurs inside recombination. 

The drag force induced on CMB photons by the scattering 
in the hydro gen recombination lines was already discussed by 
lLoebl d200ll) . They showed that due to the low population of 
the excited levels, the characteristic time over which the pe- 
culiar velocity of the gas is damped due to the drag force on 
the hydrogen atoms is much longer than the Hubble time at 
that epoch. We repeated this computation and found that it is 
at least five orders of magnitude longer. Thus the drag force 
can be safely neglected when computing the effect of coherent 
scattering in t hese lines on the power spec trum of the CMB, as 
it was done in Zaldarria ga & Loebl J2002I) and BHMS04. 

3.1. Widths of the recombination lines 



Figure [5] shows the optical depth i n the Ha line, tog ether 
with the Thomson visibility function JSunvaev & Zeldovichl 
T / j-(f7) = (dTT/drf)exp(-TT), and the function 
exp(-T7-). Here, tt is the Thomson optical depth and r\ is the 
conformal time. As we can see from this figure, the maximum 
optical depth for this transition is reached beyond the peak of 
the visibility function, around z ~ 1400, but there the Thomson 
optical depth is very large. This is why the maximum contribu- 
tion to the effect we are discussing comes from lower redshifts 
(z ~ 1000). As we show in the next section, the generation 
of new anisotropies is conected to the term tho exp(-Tr), also 
shown in Fig.|5] 



Before hydrogen is significantly recombined in the Universe, 
the optical depth for electron scattering is very high, so ev- 
ery scattering leads to a broadening of the line. For z = 1000, 
the thermal (Doppler) width of the lines is close to (Av/v),/, ss 
4x 10~ 5 . Subsequent electron scattering of these photons might 
increase this width up to (Av/v) e » 2 x 10~ 3 . The characteris- 
tic width of recombination in redshift space, (&/z) rec , is much 
bigger than this electron broadening (from Fig. [5] the visibility 
function has a width of 6z ~ 80, so (Sz/z) lec » 1/14). On the 
other hand, the optical depths in the lines under discussion are 
also changing slowly. The redshift interval (6z/z)Ha in which 
tho changes significantly (i.e. drnal^Ha ~ 1) is about 0.07 
for z = 1000. Thus, we expect that both Doppler and electron 
broadening of the lines will not change our effect significantly. 
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We have checked this point by repeating the computations of 
this paper for different widths for the lines. 

The opacity in each line (m — > n) is computed as 



• - dTL <D( \ 

TL = -T- = Trim VW, T]L) 

drj 



(3) 



where Virf, rj L ) is a profile function of area unity, and centered 
at t]i, the conformal time corresponding to the redshift that we 
are observing (1 + z a bs - v nm lv b s )- For simplicity, we have 
adopted here a normalized Gaussian profile, so 



exp(- 



2o-? 



4 



2xo* 



where 07, is the width of the Gaussian. We have explored three 
different values for o~lIt]l, namely 10~ 3 , 5 x 10~ 3 and 1(T 2 , but 
we found similar results for all of them, as expected. All values 
quoted in this paper correspond to the case (TlIi~\l = 0.005. 

One important consequence that can be determined from 
here is that every frequency observed today corresponds to a 
given redshift z. With electron scattering, we are averaging 
all the effects well before the peak of the visibility function, 
and we are losing information about the redshift in which a 
given part of the signal was produced. However, the study of 
these lines will permit us to check the velocities and optical 
depths at any stage of recombination. In order to do this, it will 
be necessary to optimise the observing frequencies of the de- 
tectors, and to have observing bandwidths narrower than the 
widths of the features (we will see below that typical widths 
are Av/v as 0.1). Unfortunately, present day and planned ex- 
periments (like Planck) have widths of the channels broader 
than this, and even broader than the width of recombination, so 
these effects are averaged inside the observing bandwidths. But 
in principle, future ground based experiments or experiments 
like CMBPol might be adapted to the width and frequency de- 
pendence of the features that we are discussing, and it could be 
possible to trace using this tomography the overall behavior of 
recombination. 

3.2. Coherent scattering in hydrogen lines during 
recombination 

In this subsection, we present the equations and the method 
of computation of the effect of coherent scattering in hydro- 
gen lines. To perform the computations for this paper, we have 
used the code from BHMS04, which was a modification of 
CMBFAST including the presence of a resonant line at a given 
redshift. We will follow here their notation. 

In the conformal Newtonian gauge, the Boltzmann equa- 
tion for the evolution of the k-mode of the photon tem- 
perat ure fluctuation can b e formally integrated to give 
( Seliak & Zaldai-riaga Hl996l) 



A T (k,T] ,fJ.) = I dr\e 
Jo 



!fc/i(j7-i7o) e -T x 



f ( A™ + ifJVb) + if>- ikfiift 



(4) 



where we have explicitly dropped the polarization term, which 
contributes at most a few percent of the total signal and will be 



neglected here. A7-0 accounts for the intrinsic fluctuations; v& 
is the velocity of baryons; <p and ijj are the scalar perturbations 
of the metric in this gauge, and ft = k ■ n. The total optical 
depth is defined here as t(t]) = tdrj, and for our problem, it 
has contributions from both Thomson and coherent scattering. 
Thus, using equation^ we have that for a given line L, the total 
optical depth is f = t j + Tl = an e {rj)o-T + TnmP(j\\ >1l)- 

Using this modification inside the CMBFAST code, and ne- 
glecting the drag force induced by this process on the CMB 
photons, it is possible to compute the effect of the coherent 
scattering on the power spectrum. The important point for us is 
that, as shown by HMS05, the effect of lines on the CMB spec- 
trum is amplified due to correlation with the existing radiation 
field. The presence of coherent scattering in lines produces an 
angular redistribution of the CMB photons. For the lines con- 
sidered in BHMS04, the net effect is generation of anisotropics 
at large scales due to Doppler motion of scattering atoms, and 
suppression of power at small scales. Summarizing, the ampli- 
fication effect consists in that the change in the power spectrum 
{5C( = C'p e - C() due to a given transition m — > n is propor- 
tional to r„„, (and not to r 2 ,,,) for small values of the optical 
depths. Thus, in the optically thin limit we have 



6C e = r nm ■ Ci(t) + tL ■ C 2 {() + 0(tL) 



(5) 



were the equation for the functions C\(t) and C2(f) is given in 
Appendix A of BHMS04 (note that in that equation, the po- 
larization terms were neglected). In the limit of high-ff values, 
they obtained the well-known result of C\{() — > -2Q. 

Eq. [5] can be used in our problem, in which the scatter- 
ing line is embedded inside recombination, although there are 
some differences to the case in which the line is placed between 
recombination and us. As was shown in BHMS04, the genera- 
tion and blurring terms tend to cancel each other at a multipole 
value I ~ (770 - J] cs )/(j] cs - j] rec ), with T] cs , j] dec and 770 the con- 
formal times at the coherent scattering, at recombination and 
the present day, respectively. From this, one can expect that 
if coherent scattering takes place not far from recombination, 
both generation and blurring will cancel over a wider range of 
multipoles. We illustrate this in Fig. [6] where we show the de- 
pendence of the linear term (C\{()) for four different cases of 
a hypothetical line placed at redshifts z = 850,950, 1050 and 
1 150. Thus, in other words, the presence of this strong decrease 
in the power at angular scales t % 200 is direct evidence that the 
line is located at the epoch of recombination in the Universe. 
Unfortunately, because of this decrease, the detection at large 
angular scales of this signal will be more difficult. 

Only at redshifts lower than the redshift of recombination 
is the blurring term dominant at high multipoles: when the line 
is well within the peak of the Thomson visibility function, both 
blurring and generation terms cancel each other, at least to a 
level of a few percent. This is seen in Fig. [6] so for redshifts 
in the tail of the visibility function (z = 850) the linear term 
is equal to -2Q at high multipoles. But if we move to higher 
redshifts, this behavior disappears. 

For illustrational purposes, we discuss how the visibility 
function is modified when we consider coherent scattering in 
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Fig. 6. Linear term of the correction to the angular power spec- 
trum originated from coherent scattering in a line located at the 
specified redshift. These curves have to be multiplied by the 
optical depth in the line to obtain the first order correction to 
6C(. Thick lines correspond to positive values, and thin lines 
to negative ones. For comparison, it is also shown the primary 
CMB power spectrum (dotted lines). Note that at high-f, the 
correction is roughly -2Q for the case z = 850 or z = 950, as 
we would expect for a line located between us and recombina- 
tion. This property disappears when we go to higher redshifts, 
because in that case the line is located inside recombination 
where radiative viscosity and thermal conductivity are impor- 
tant. [See movie in the electronic version]. 

the hydrogen lines. In the optically thin limit, we can write to 
first order in the optical depth (t„,„) that 

'Vfoifc) = te~ T * [l-r,„„y[(7,)]^ r +r„ ffl ;P(7,) e - r r + 0(TL)(6) 
where we have defined the area function of the profile, 3ttr\), 



f 

J 7] 



dim) 



From here, we see that to first order in the optical depth, we 
have two terms giving us the correction to the Thomson visi- 
bility function: one associated with suppression 
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Fig. 7. Correction to the visibility function in the presence of 
coherent scattering in the Ha line. To first order in the optical 
depth of the line (r nm ), there are two terms. One is associated 
with suppression (eq.0, and the other one with the generation 
of new anisotropies (eq. [8J. Both terms are plotted (dot-dashed 
and dashed lines, respectively) together with their sum (solid 
line), for the case of observing the Ha line placed at z = 1000. 
For comparison, we also show the Thomson visibility function 
(dotted line). 



and another connected to generation of anisotropies: 

<P(n)e~ TT 



(8) 



% w = -T nm 3[{rj)'V T = -T nm J[{ri)f T e Tj 



(7) 



^ 'gen — T\ 

In Fig. we display these two terms for a case in which we 
are observing the Ha line placed at z = 1000 (i.e. we are ob- 
serving at vo = 456 GHz), and with an instrumental width 
of o~l/t]l = 0.005 . The dot-dashed line corresponds to the 
suppression term, whereas the dashed line shows the genera- 
tion one. The former is proportional to the Thomson visibility 
function for redshifts larger than the redshift being probed by 
the line, since only anisotropies located behind the line can be 
blurred. On the other hand, the generation of anisotropies is lo- 
cated at the position of the line, and its amplitude is weighted 
by the amount of free electrons situated between the line and 
the observer (factor exp(-Tj-) in eq.[S}. The curve T Ha exp(-r r ) 
was shown in Fig. [5] 

3.3. Results for the considered lines 

Our main result for the nine transitions under discussion can be 
summarized in Table|2] where we present the angular scale and 
the redshift at which we obtain the maximum signal for each 
of the lines. From these values, we see that the best angular 
scale to look for this effect is placed in the region of the third 
acoustic peak, and that the maximum signal comes from red- 
shifts z ~ 1010 - 1060. These values lie close to the peak of the 
visibility function, but in the low-redshift tail. We also present 
the full widths of the regions around the maxima where the sig- 
nal is within 70% of the peak value, both in multipole (AC) and 
frequency (Avo). We note that Avq/vq is of the order of 10%, 
which is hundreds of times larger than the Doppler broadening 
of the lines. This will help in the detection of these features. 



J. A. Rubino-Martm: Cosmological hydrogen recombination lines and CMB 



9 



Table 2. Maximum signal (in absolute value) obtained for each of the lines within the considered cosmological model. It is shown 
the angular multipole, the redshift and the observed frequency (today) at which we have maximum effect. We also present the 
(full) widths At and Av ( > around the peak value, defined as the regions where the signal is greater than 70% of the maximum 
value. The last two columns show the amplitude of the effect in temperature QiK) and the relative value of the temperature to the 
primordial CMB power spectrum. All these maximum deviations are negative (i.e. the correction to the power spectrum SCt at 
those maxima is negative). 
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We also show that the width of the feature in multipole space 
is of the order of 140, so a full sky coverage is not necessary to 
detect the signal. 

The last two columns in Table [2] show the amplitude of 
the signals at these maxima, both in temperature and the frac- 
tional change in temperature with respect to the intrinsic CMB 
power spectrum. Although the values for -yVQI/Q are small, 
the present day sensitivity for the detection of angular fluctua- 
tions is better than that for spectral measurements, so in prin- 
ciple these features should be easier to detect than the distor- 
tions of the spectrum. We now investigate in more detail the 
angular and frequency dependences of this effect. We will con- 
centrate here on the Ha line, which gives the strongest signal. 
Predictions for the other transitions can be derived in a similar 
way. 

Fig. shows the angular dependence of the effect, for the 
case of Ha line, and for redshift values z = 900, 1010 and 
1150, which correspond to (redshifted) present day frequen- 
cies of 507, 452 and 397 GHz, respectively. Given the small 
values of the optical depth, the linear term gives the dominant 
contribution to the effect, so the shape of the plot is similar 
to Fig. |6] but the amplitude is modulated by the optical depth 
for the transition at every redshift. As shown in Table |2] the 
maximum signal is obtained at redshift z = 1010. Prior to this 
redshift, the shape of the 6Cc mimics that of the primordial 
C( at high multipoles, because the suppression term dominates 
here (and hence the 6C( are negative). When the line is placed 
at higher redshifts, well inside the recombination, the cancel- 
lation of the generation and blurring terms at high multipoles 
produces a number of regions where the 5C{ quantities have 
different signs. The exact position of the zeros strongly depends 
on the particular cosmological model, and as pointed above, it 
could be used to investigate in detail the shape of the transfer 
function at all redshifts. 

In Fig. |5] we present the same redshift slices as in Fig. [6] 
but now showing \6Cc\/Cc for the Ho- transition. The relative 
value of the distortion of the power spectrum is of the order of 
10~ 5 - 10~ 4 in band power (~ 10~ 2 in temperature). It should be 
also noted that \SCe\/C( is growing at large multipoles, and ap- 



proaches the asymptotic value of -2thu (shown here as a dot- 
ted line), although the multipole region where we reach this be- 
havior depends on the considered redshift slice. When the line 
is placed at higher redshifts inside recombination, the asymp- 
totic behavior is reached at higher multipoles in the damping 
tail region of the spectrum. Moreover, when the redshift is high 
enough, this asymptotic region disappears, and we find a com- 
plicated structure of positive/negative features due to the afore- 
mentioned cancellation. Thus, this method co uld give us infor- 
mation about the 'Silk damping' JSilkLll968l) . or processes of 
dissipation of short wavelength adiabatic perturbations due to 
radiative viscosity and thermal conductivity during the recom- 
bination. 

We now consider the frequency dependence of the effect. 
Figure^|shows the relative change in the power spectrum due 
to the presence of the hydrogen recombination lines of Balmer 
series, as a function of the present day redshifted frequency. For 
each Balmer line, we show three different panels correspond- 
ing to three different angular scales (I = 20, 300 and 800). As 
pointed above, the signal has a very peculiar frequency depen- 
dence, and for example at angular scales of I — 20 and t = 800, 
we can see a characteristic change of sign of the effect, related 
to the positions of the cancellations of the generation and blur- 
ring terms. For the most intense line, Ha, we show five more 
angular scales in Fig. ^2 where we explicitly use a linear scale 
in the vertical axis to show the change of sign. The case of 
t = 873, which is giving the maximum signal (in absolute 
value) is shown as a solid line in that figure. The frequency de- 
pendence of the effect is unique and completely different from 
that of foregrounds (e.g. dust emission), and thus it might be 
used to separate these signals from other contaminants, as we 
will see below. 

The figures for the Paschen and Brackett lines are similar to 
those for Balmer lines. For illustration, in Figure^]we present 
the angular and frequency dependence for the Ba and the Pa 
lines. The amplitude of the signal for all the transitions consid- 
ered in this paper is presented in Tabled 

Summarizing this subsection, we conclude that the best 
strategy to observe these features and perform this tomography 
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Fig. 8. Correction to the power spectrum due to coherent scat- 
tering in the Ha line. We plot three curves for three different 
redshifts, namely z = 900, 1 150 (which correspond to 507 GHz 
and 397 GHz, respectively), and also z = 1010 (452 GHz), 
which is the redshift (frequency) at which we have the max- 
imum amplitude for signal, at multipole t = 873. Note how 
the z = 900 case mimics the shape of the power spectrum at 
high multipoles, but with negative sign, and how this behavior 
disappears as we go to higher redshifts. Observing at differ- 
ent frequencies with sufficiently narrow spectral bands, we can 
perform tomography of the recombination epoch using the Ha 
line. 

would be to use bandwidths (Av/v) of the order of 10%, al- 
though bandwidths of 1 % would trace this signal much better. 
These values are narrower than the width of the visibility func- 
tion in the line, but still wide enough to keep enough photons. 
It is clear that the major contribution to the effect during the 
epoch of recombination is coming from high I multipoles, and 
that the widths of the features is of the order of A£ ~ 140, so 
this opens the possibility to use not only satellites but balloons 
or ground based experiments to look for this signal. 

4. Contamination from other effects 

In this section we discuss how the presence of foregrounds 
and the effect associated with Rayleigh scattering might affect 
the study of the lines. We shall adopt the middle-of-the-road 



Fig. 9. Relative correction QSCe |/ Q) to the power spectrum due 
to scattering in the Ha line, as a function of angular scale {€). 
Thick lines correspond to positive values whereas thin lines 
represent negative ones. The distortions are larger at small an- 
gular scales, approaching the asymptotics of -2t Ho (dotted 
line) for redshifts before the peak of the visibility function. The 
region where this asymptotic behavior is valid changes with 
redshift, and for redshifts well inside the recombination, it dis- 
appears. This behavior could give us information about the pro- 
cess of dissipation of short wavelength adiabatic perturbations 
due to radiative viscosity and thermal conductivity during re- 
combination. 

model of iTegmark et alJ (2000, hereafter TOO), and we con- 
sider the contribution of five foreground sources, namely syn- 
chrotron radiation, free-free emission, dust emission, thermal 
SZ effect associated with filaments and superclusters of galax- 
ies, and Rayleigh scattering. Point sources are not considered 
here, so we assume that resolved sources can be extracted from 
the maps, and that the contribution of unresolved sources can 
be lowered down to roughly the noise level of the observing 
instrument. 

For the first three components, synchrotron, free-free and 
dust, the angular dependence of the power spectra was approx- 
imated by a power law. The frequency dependence for free-free 
and synchrotron was taken to be a power law, whereas a modi- 
fied black-body spectrum was used for the dust component. For 
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Fig. 12. Angular power spectrum arising from the Pa (top row) and Ba (bottom row) hydrogen lines during recombination, as a 
function of the redshifted frequency (left column) and the angular multipole (right column). The solid lines in all panels refer to 
the cases in which the signal is maximum, according to Table[2] 



the SZ effect, the frequency dependence is well-known. Details 
of the modeling of the angular dependence can be found in TOO. 
This neglects the contribution from the SZ effect generated by 
resolved clusters of galaxies, which we assume can be removed 
from the maps. For the Rayleigh scattering effect modeling, 
we follow BHMS04. As discussed above, the scattering cross- 
section of this process has a very strong dependence on the 
frequency, 



<tr(v) = o~tv 



z 

V k>2 



fu 



(v 2 lk -v 2 )) 





1 \ 




V 


oc 




/ 


V\2 



where the last proportionality holds for v <c v\2- Thus, the 
net effect on the power spectrum will also show this depen- 
dence (6Cf oc tr oc v 4 ). Unfortunately, this frequency depen- 
dence mimics the spectrum of dust emission from local dust 
and bright extragalactic star-forming galaxies, so the separa- 
tion of these effects could be difficult. 

Fig- El shows the expected level, prior to any removal, of 
the discussed foregrounds when looking for the maximum of 
the Ha line. We considered here the case of an ideal experi- 
ment measuring at 430 GHz and 130 GHz, and with an instru- 
mental bandwidth of Av/v = 10~ 2 , which is comparable to the 



electron width of the line under discussion. It is clear that the 
main diffuse contamination will come from dust emission, but 
there are other two signals which will be larger than the Ha 
contribution, namely the thermal SZ emission from filaments 
and the Rayleigh scattering. For the first one, the signal has a 
completely different frequency and angular behavior from the 
one we are considering, so these peculiarities could be used to 
separate both components. 

For the case of Rayleigh scattering, the frequency behavior 
is also completely different from the case of resonant scatter- 
ing. The angular dependence has some similarities in the high 
multipole range (damping tail) if we probe a line in the low 
redshift tail of recombination, because here both effects are 
proportional to -2tQ . However, as we probe higher redshifts, 
the behaviors become different also in this domain. In addi- 
tion to this, there is another difference between the Rayleigh 
scattering and the coherent scattering in lines, which is related 
to the different evolution of the optical depth of each process 
during recombination, and which is superimposed on the pre- 
vious effect. To illustrate it, we show in Fig. [21 the redshift 
dependence of the population of electrons, neutral hydrogen 
atoms and the population of hydrogen atoms in the level 2p, 
all computed using our code. These three variables (n e , tim 
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Fig. 10. Angular power spectrum arising from the Balmer re- 
combination lines of hydrogen during cosmological recombi- 
nation, as a function of the redshifted (observed) frequency. We 
plot the values at t — 20 (top panel), at t — 300 (middle panel), 
and at t — 800 (bottom panel). Thick lines corresponds to re- 
gion where the 5C[ are positive, whereas thin lines correspond 
to negative values. 



and nip) are proportional to the differential optical depth t for 
the Thomson, Rayleigh and coherent scatterings, respectively. 
From this, it is clear that the redshift interval dz/z in which the 
optical depth for Rayleigh scattering changes significantly is 
much larger than the corresponding one of the Ha line scatter- 
ing. Thus, instruments with broad bandwidths will dilute the 
signal from the lines, but add up the signal from Rayleigh scat- 
tering. This also could be used to separate them, because using 
a broad-band instrument we can isolate the Rayleigh scattering 
component, and then subtract it. 

One interesting physical remark is that the evolution of 
the population of the 2p level in Fig. El is practically pro- 
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Fig. 11. Same as Figure flOl but only for the Ha line which gives 
the maximum contribution. The vertical axis is now presented 
in a linear scale so we can explicitly see the change of sign. We 
plot values for five multipole scales, t = 100,350,750, 1200, 
and also I = 873, which is the multipole at which we have the 
maximum signal at v = 452 GHz. [See movie in the electronic 
version] . 

portional to the electron density in this redshift range. This 
is what we expected, as we can immediately see from Eq. |2 
(ri2p °c n 2 e expl^/feT^)}) and from the evolution of the elec- 
tron density (x e oc (1/z) expf-^/feTe)}, see ZKS68). 

5. Discussion and conclusions 

We have studied here the imprint of cosmological hydrogen re- 
combination lines on the power spectrum of the CMB. To this 
end, we have developed a code that follows the evolution of the 
population of the levels of the hydrogen atom, separately treat- 
ing those levels with different angular momentum. From here, 
we have obtained the optical depths associated with coherent 
scattering in these lines during the epoch of recombination, and 
these numbers were used to compute the effect on the angular 
power spectrum of the CMB. 

These changes are small (0. l/jK-03fiK), but could be sepa- 
rated from other effects due to their peculiar frequency and an- 
gular dependence. Unfortunately, the line giving the maximum 
signal (the Ha line) is placed in a frequency domain where the 
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Fig. 13. (a) Maximum contribution of the coherent scattering in the Ho- line to the angular CMB power spectrum, as seen by an 
hypothetical experiment measuring at 430 GHz and 130 GHz. The upper thick solid line gives the reference power spectrum, (b) 
Contribution of foregrounds at these frequencies. All thin lines refer to the foreground model as quoted in TOO: free-free emission 
is shown as dotted line; synchrotron emission as solid line; the thermal SZ effect produced in filaments as dot-dashed line and the 
thermal dust emission is shown by dashed lines. Rayleigh scattering contribution is shown as a thick long-dashed line, whereas 
the cosmic variance contribution associated to this effect is shown as a thin long-dashed line. 



contribution from dust emission, the tSZ effect produced in fil- 
aments and Rayleigh scattering are important, so it will be nec- 
essary either to look for regions with low contamination (for 
the case of dust), or to use component separation techniques 
to reach this signal. The important point here is that the sig- 
nal under discussion has a very characteristic behavior, both in 
frequency and in angular scale, which is completely different 
from any of the above contaminants. 

One of the most important properties of this signal is that 
each observing frequency is associated with one redshift, so ob- 
servations of these signals at different frequencies might give 
information about the amplitude of the fluctuations and cor- 
responding peculiar velocities of matter at different redshifts 
during the epoch of recombination. 

In this pap er, the best-fit cosm ological model to the WMAP 
lst-year data (Bennett et al., 2003) was adopted to make pre- 
dictions about the shape and the intensity of the signal we ex- 
pect to measure. With the announced sensitivities of future ex- 
periments, these signals should be detected, and thus we could 
try to infer some information on the values of the cosmologi- 
cal parameters. The inferred constraints from the detection of 
these lines will be clearly independent of those obtained using 
the model-fitting to the power spectrum of the CMB. This ad- 
ditional information could be used to break some degeneracies 
in the parameters in two ways. First, as pointed out above, ob- 
servations in narrow spectral bands of these signals can give us 



information about specific redshift slices during the process of 
recombination in the Universe. The shape and positions of the 
peaks and zeros of the signal presented in Figs.[S]and[9]reveal 
the details of the baryon transfer function at each considered 
epoch, so they could be used as an additional test of the cosmo- 
logical model. For a given experiment with a given bandwidth, 
it is straight-forward to derive the prediction for 6C(, so model 
fitting could be used here, in the same way as it is used today 
to extract parameters from the power spectrum. Secondly, the 
determination of the amplitude of this effect at a given redshift 
could also be used to set an independent constraint on Q^/i 2 . 

Using this signal, future narrow-band spectral observations 
might permit one to study all the details of the transfer function 
describing the evolution and viscous damping of perturbations 
in the epoch of recombination. 
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